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We study the critical dynamics of the generalized Smoluchowski-Poisson system (for self- 
gravitating Langevin particles) or generalized Keller-Segel model (for the chemotaxis of bacterial 
populations). These models [Chavanis & Sire, PRE, 69, 016116 (2004)] are based on generalized 
stochastic processes leading to the Tsallis statistics. The equilibrium states correspond to poly- 
tropic configurations with index n similar to polytropic stars in astrophysics. At the critical index 
nz — d/(d — 2) (where d > 2 is the dimension of space), there exists a critical temperature Q c (for 
a given mass) or a critical mass M c (for a given temperature). For O > O c or M < M c the system 
tends to an incomplete polytrope confined by the box (in a bounded domain) or evaporates (in an 
unbounded domain). For O < C or M > M c the system collapses and forms, in a finite time, a 
Dirac peak containing a finite fraction M c of the total mass surrounded by a halo. We study these 
regimes numerically and, when possible, analytically by looking for self-similar or pseudo self-similar 
solutions. This study extends the critical dynamics of the ordinary Smoluchowski-Poisson system 
and Keller-Segel model in d = 2 corresponding to isothermal configurations with ns — > +oo. We 
also stress the analogy between the limiting mass of white dwarf stars (Chandrasekhar's limit) and 
the critical mass of bacterial populations in the generalized Keller-Segel model of chemotaxis. 



I. INTRODUCTION 



For a long time, statistical mechanics was restricted 
to systems interacting via short-range forces. For exam- 
ple, the case of self-gravitating systems is almost never 
considered in standard textbooks of statistical mechan- 
ics and these systems have been studied exclusively in 
the context of astrophysics. In the sixties, Antonov 
@, Lynden-Bell @ and Thirring @ realized that self- 
gravitating systems have a very special thermodynamics 
marked by the non-equivalence of statistical ensembles 
(microcanonical, canonical, grand canonical,...). This is 
related to the non-additivity of the energy and to the 
presence of negative specific heats in the microcanoni- 
cal ensemble. Furthermore, these systems experience a 
rich diversity of phase transitions (microcanonical and 
canonical first order phase transitions, zeroth order phase 
transitions,...) associated with their natural tendency to 
undergo gravitational collapse 0, M. Recently, several 
researchers have started to consider the dynamics and 
thermodynamics of systems with long-range interactions 
at a more general level (see the books || [7J and refer- 
ences therein) and to discuss the numerous analogies (and 
differences) between these systems: self-gravitating sys- 
tems, two-dimensional vortices, neutral and non-neutral 
plasmas, the HMF model, free electron lasers, Bose- 
Einstcin condensates, atomic clusters, chemotaxis of bac- 
terial populations etc. These analogies have also sug- 
gested interesting experiments. For example, in the 
physics of ultra-cold gases, some authors H have pro- 
posed to generate an attractive 1/r interaction between 
atoms by using a clever configuration of laser beams. 
This leads to the fascinating possibility of reproducing, in 
the laboratory, the isothermal collapse (in the canonical 
ensemble) of a self-gravitating Fermi gas M no) leading 



to a "white dwarf star" . These examples illustrate the 
importance of studying the statistical mechanics of sys- 
tems with long-range interactions at a general level and 
to develop the analogies between different systems that 
may seem a priori of a very different nature. 

In a scries of papers @ g^, [y| [u], |l|, ||, , we have 
investigated the dynamics and thermodynamics of a sys- 
tem of self-gravitating random walkers. The basic idea is 
to couple the usual Brownian motion (as introduced by 
Einstein and Smoluchowski) to the gravitational interac- 
tion. In our general model |L7fl , the microscopic dynam- 
ics of the particles is described by N coupled stochastic 
equations including a friction force and a stochastic force 
in addition to the gravitational interaction. The fric- 
tion force and the stochastic force model the interaction 
of the system with a thermal bath of non-gravitational 
origin. Then, the proper statistical description of this 
dissipative system is the canonical ensemble. In order to 
simplify the problem, we have considered a strong fric- 
tion limit in which the motion of the particles is over- 
damped. We have also considered a mean field approxi- 
mation which becomes exact in a proper thermodynamic 
limit N — > +oo in such a way that the volume V ~ 1 is 
of order unity and the coupling constant G ~ 1/N goes 
to zero (alternatively, we can consider that the mass of 
the individual particles scales like m ~ 1/N so that the 
total mass M ~ Nm and the gravity constant G re- 
main of order unity). These approximations lead to the 
Smoluchowski-Poisson (SP) system. The steady states 
correspond to isothermal distributions associated with 
the Boltzmann statistics. When coupled to the Poisson 
equation, we obtain density profiles similar to isother- 
mal stars in astrophysics |18| , |l9f . In the course of our 
study, we realized that the SP system is isomorphic to 
the standard Keller-Segel (KS) model (20| |I[ introduced 



Index 



Temperature 



Bounded domain 



Unbounded domain 



• Evaporation JL6| : 

asymptotically free normal 

diffusion (gravity negligible) 

• Collapse: 

pre-collapse and post-collapse as 

in a bounded domain |ll|, |l2], Pq] 



T> T c 



Metastable equilibrium state 

(local minimum of free energy) : 

box-confined isothermal sphere [11, 12, 35| 



T <T C 



Self-similar collapse with a — 2 |11, 12, 39] 

followed by a self-similar post-collapse leading 

to the formation of a Dirac peak of mass M Jl3| 



e>e c 



< n < m 



Equilibrium state: 
box-confined (incomplete) polytrope |5l| 



e<e c 



Equilibrium state: 
complete polytrope (compact support) felj 



Equilibrium state: 

complete polytrope 

(compact support) fell] 



e>e c 



n-j < n < oo 



Metastable equilibrium state 

(local minimum of free energy) : 

box-confined polytropic sphere |5l| 



e<e c 



Self-similar collapse with a — 2n/(n — 1) [ pl| 
followed by a post-collapse leading to the 
formation of a Dirac peak of mass M [N] 



• Evaporation [P]: 

asymptotically free anomalous 

diffusion (gravity negligible) 

• Collapse: 
pre-collapse and post-collapse 
as in a bounded domain felj 



e>e c 



n = n 3 



Equilibrium state: 
box-confined (incomplete) polytrope |5l| 



Self-similar evaporation 
modified by self-gravity [P] 



e<e c 



Pseudo self-similar collapse 

leading to a Dirac peak of 

mass (e/9 c ) d/2 M + halo [P]. 

This is followed by a post-collapse 

leading to a Dirac peak of mass M [N] 



Collapse [N] 



e = e c 



Infinite family of steady states [P] 



Infinite family of steady states [P] 



TABLE I: Summary of the different regimes of the GSP system in d > 2 with references to the physical literature ([P]: present 
paper; [N]: not done). The case of negative indices is considered in | p2| . The links to the mathematical literature are indicated 
in the main text. Note: for (n = oo, T > T c ) and for (713 < n < 00, O > O c ) in a bounded domain, the system can either reach 
a metastable equilibrium state or collapse depending on a notion of basin of attraction (see |ll[ for more details). 



in mathematical biology to describe the chemotaxis of 
bacterial populations (22j . The SP system and the KS 
model have now been extensively studied by physicists 
0, |1 p|, P, p|, pi, |j~7l, |23| a nd applied mathcmati- 



30l |3U m, HI ||, ||, 1 



10|, y,]42f|4|r^ri|r^,^ with different 
methods and motivations. 

We have also studied a generalized Smoluchowski- 
Poisson (GSP) system [see Eqs. |L3])-([l4|) of this pa- 
per] including an arbitrary barotropic equation of state 
P(p). This model has been introduced by Chavanis in 
p8| . The GSP system can be viewed as a generalized 
mean field Fokkcr-Planck equation (for a review of non- 
linear Fokker-Planck equations, see M, pG|). It can be 
obtained from generalized stochastic processes and it is 
associated with a notion of effective generalized thermo- 
dynamics (E.G.T). These equations can also provide a 
generalized Keller-Segel (GKS) model of chemotaxis with 
a density dependent diffusion coefficient |50[ . For an 
isothermal equation of state P = pksT /m, we recover 
the standard SP system and KS model (with appropri- 
ate notations). Apart from the isothermal equation of 



state, the GSP system and GKS model have been stud- 
ied for: (i) a polytropic equation of state P = Kp 1 |5l]] 
(ii) a logotropic equation of state P — A In p [J52| ( iii) a 
Fermi-Dirac equation of state P = Pf.d.(p) H53|, |54| (iv) 
an equation of state P = — Tp max ln(l — p/p ma ,x) tak- 
ing into account excluded volume effects p5[ . These are 
standard equations of state introduced in astrophysics 
and statistical mechanics so that it is natural to consider 
these equations of state in connexion to the GSP system 
and GKS model. 

Specializing on the polytropic equation of state P = 
Kp 1 with 7 = 1 + 1/n pa, the steady states of the 
GSP system correspond to polytropic distributions asso- 
ciated with the Tsallis statistics fjq]. When coupled to 
the Poisson equation, we obtain density profiles similar to 
polytropic stars in astrophysics JLq , [L9j. For d > 2, there 
exists a critical index 74/3 = 2(d—l)/d, i.e. n^ = d/(d—2) 
pl[ . For < n < 713, the GSP system relaxes towards 
a stable steady state with a compact support, similar to 
a classical white dwarf star (classical white dwarf stars 
are equivalent to polytropcs with index n = 3/2 in d = 3 
p7|). For n > 713, there is no stable equilibrium in an 



Index 


Temperature 


Bounded domain 


Unbounded domain 


n = oo 


T>T C 


Equilibrium state: 
analytical solution |12| 


Self-similar evaporation 
modified by self-gravity |lq] 


T <T C 


Pseudo self-similar collapse 

leading to a Dirac peak of 

mass (T/T C )M + halo Jl|, 0. 

This is followed by a post-collapse 

leading to a Dirac peak of mass M [N] 


Collapse [N] 


T = T C 


Self-similar collapse leading to 

a Dirac peak of mass M with 

exponential growth of p(0, i) |12l |47|] 


Self-similar collapse leading to 

a Dirac peak of mass M with 

logarithmic growth of p(0,t) [fL SI] 


< n < oo 


T>T C 


Equilibrium state: 
box-confined (incomplete) polytrope fplf 


Equilibrium state: 
complete polytrope (compact support) 


T<T C 


Equilibrium state: 
complete polytrope (compact support) |5l|] 



TABLE II: Summary of the different regimes of the GSP system in d — 2. In d ■■ 
statistical equilibrium state so that there is no evaporation or collapse ]12|, B3, p 



1, the GSP system always relaxes towards a 



unbounded domain so that the system can either collapse 
or evaporate (see Fig. [13] for an illustration). These dif- 
ferent regimes have been studied in plj |. For n — 713, the 
dynamics is critical. At this index, there exists a criti- 
cal mass M c (d) (for a given polytropic constant K) Q 
which is connected to the Chandrasckhar mass of rcla- 
tivistic white dwarf stars (ultra-relativistic white dwarf 
stars are equivalent to polytropes with index n = 3 in 
d = 3 [p8| ). The object of the present paper is to study 
numerically and, when possible, analytically this criti- 
cal dynamics. For M < M c , we find that the system 
evaporates and we construct a self-similar solution. For 
M > M c , we find that the system collapses. In a finite 
time ton, it forms a Dirac peak with mass M c surrounded 
by a halo that has a pseudo self-similar evolution. For 
d = 2, the critical index 713 — ► +00 so that we recover 
the case of isothermal spheres whose dynamics is known 
to be critical in d = 2 Q . 

When we apply this model in the context of chemo- 
taxis p9|, we find the existence of a critical mass M c (d) 
at the critical index n 3 = d/(d—2). For d = 2, we recover 
the well-known result M c (d = 2) = 8tt obtained within 
the standard Keller-Segel model (see [p(| and references 
therein) and for d = 3, the critical mass associated with 
the GKS model is M c (d = 3) = 202.8956... (in usual 
dimensionless variables). This is similar to the Chan- 
drasckhar limiting mass of white dwarf stars. The exis- 
tence of a limiting mass for bacterial populations at the 
critical index 713 and its connexion to the Chandrasckhar 
mass was pointed out in [[39[ ^l) (and implicitly in plf ). 
This is another illustration of the numerous analogies 
that exist between self-gravitating systems and bacterial 
populations p3[ . 

The paper is organized as follows. In Sec. |l[ we 
briefly recall the connexion between white dwarf stars 



III 



we recall the ba- 



and gaseous polytropes. In Sec. 
sic properties of the SP and GSP systems and describe 
the behavior of the solutions depending on the index n 
and the dimension of space d. As the problem is very 
rich, involving many different cases (~ 30), a summary 
of previously obtained results, completed by new results 
and new discussion, is required to understand the place 
of the present study in the general problem (see also Ta- 
bles I and U for an overview). Then, we consider more 
specifically the particular index n = n-j which presents 
a critical dynamics that was mentioned, but not stud- 
ied, in our previous paper pl[. 



IV, we show 



m our previous paper [bi|. In Sec 
that this critical value can be understood from a sim 
pie dimensional analysis. In Sec. M, we study the crit- 
ical collapse dynamics and extend the results obtained 
in d = 2 for isothermal (n = +00) systems |I5] to the 
case of critical polytropes (n = 71.3) in d > 2. In Sec. VI, 
we study the evaporation dynamics in unbounded space. 
We show that for n > 713, self- gravity becomes negligible 
for large times so that the evaporation is eventually con- 
trolled by the pure (anomalous) diffusion. For n = 713, 
gravity remains relevant at any time so that there exists 
a self-similar solution for which all the ter ms o f the GSP 
system scale the same way. Finally, in Sec. VII, we trans- 



pose our main results to the context of chemotaxis using 
notations and parameters adapted to this problem (this 
is to facilitate the comparison with the results obtained 
in mathematical biology). 

Our numerical and analytical study was conducted in 
parallel to a mathematical work by Blanchet et al. I64] 
who obtained rigorous results for the critical dynamics 
of the GSP system and GKS model introduced in our 
paper pi] . These two independent studies have different 
motivations and use very different methods so they are 
complementary to each other. 



II. WHITE DWARF STARS AND POLYTROPES 

In this section, we briefly recall the connexion be- 
tween the maximum mass of white dwarf stars (Chan- 
drasekhar's mass |J58|1) and the theory of self-gravitating 
polytropic spheres |18[ J19II ■ 

In simplest models of stellar structure, a white dwarf 
star can be viewed as a degenerate gas sphere in hydro- 
static equilibrium. The pressure is entirely due to the 
quantum pressure of the electrons (resulting from Pauli's 
exclusion principle for fermions) while the density of the 
star is dominated by the mass of the protons. The con- 
dition of hydrostatic equilibrium coupled to the Poisson 
equation reads 
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and the equation of state of a degenerate gas of relativis- 
tic fermions at T = can be written parametrically as 
follows [p3| 



P = A 2 f(x), p = Bx 3 



where 



An 



„.4 5 
7TTO C 



B = 



&irm i c J 'pH 
3p 



f(x) = x(2x 2 - 3)(1 + x 2 ) 1/2 + 3 sinh" 



(2) 



(3) 



(4) 



where m is the mass of the electrons, H is the mass of 
the protons and p, is the molecular weight. The function 
f(x) has the asymptotic behaviors f(x) ~ (8/5)cc 5 for 
x <C 1 and f(x) ~ 2x 4 for x 3> 1. The classical limit 
corresponds to i < 1 and the ultra-relativistic limit to 
i»l. In these limits, the white dwarf star is equivalent 
to a polytropic gas sphere with an equation of state P = 
Kp 1 . The index n of the polytrope is defined by 7 = 
I + I/7;,. In d = 3 dimensions, polytropes are self-confined 
for n < 5 and they are stable (with respect to the Eulcr- 
Poisson system) for n < 3 (for n = 3 they are marginally 
stable) . The mass-radius relation is given by @ : 



M (n-l)/np(S-n)/n _ K ( 1 + n ) (n-l)/n 

G(47r)V» W " 



(5) 



where ui n is a constant (depending only on the index n 
of the polytrope) that can be expressed in terms of the 
solution of the Lane-Emden equation [Q . 

In the classical case i< 1, the equation of state takes 
the form 



P = K lP 



5/3 



with 



A',=i|^ 



5V87T 



2/3 



i{pHf/i 



(6) 



(7) 



FIG. 1: Mass-radius relation for relativistic white dwarf stars 
at T — [p8|. The radius vanishes for a limiting mass 
Mchandra corresponding to the ultra-relativistic limit (R). 
The dashed line corresponds to the classical limit (C). 
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FIG. 2: Mass versus central density for relativistic white 
dwarf stars at T — 0. Equilibrium states only exist for 
M < M Chandra- For M — M C handra, the density profile 
is a Dirac peak. For M > Mchandra, the system is expected 
to collapse and form a neutron star or a black hole. The cor- 
responding density profiles are represented in Fig. 4 of J63]. 



Therefore a classical white dwarf star is equivalent to a 
polytrope of index n = 3/2. The mass-radius relation is 
given by 



2/3 



/r 



2V327r 2 y mG(pH) 5 / a 3 / 2 ' 



(8) 



with w 3/2 = 132.3843.... It exhibits the familiar MR 3 ~ 
1 scaling. 

In the ultra-relativistic limit x 3> 1, the equation of 
state takes the form 



P = K 2P i ' 3 , 



(9) 



with 



K 2 = 



1/3 



he 



(/X#) 4 /3 



(10) 



Therefore, an ultra-relativistic white dwarf star is equiv- 
alent to a polytropc of index n = 3. For this index, the 
relation (||) leads to a unique value of the mass 



M, 



32tt s 



1/2 fhc\ 3/2 



( M ff)S 



(11) 



with 0J3 = 2.01824.... This is the Chandrasekhar mass 



M c = 0.196701... — 



hc\ 3/2 1 



g; (f,Hy 



5.76M Q /p 2 . (12) 



Considering the general mass-radius relation of partially 
relativistic white dwarf stars (see Fig. [I]), we note that, 
for this limiting value, the radius R of the configuration 
tends to zero. This leads to a Dirac peak with mass M c . 
Thus, the Chandrasekhar mass represents the maximum 
mass of white dwarf stars (see Fig. 0) . There is no hy- 
drostatic equilibrium configuration for M > M c . 

If we extend Chandrasekhar's theory to a d- 
dimensional universe |6l| , we find that white dwarf stars 
become unstable in a universe with d > 4 dimensions (in 
d = 4, classical white dwarf stars exist for a unique value 
of the mass M = M c = 0.0U3958...h 4 /(m 2 G 2 p 3 H 3 ) and 
they are marginally stable). Therefore, the dimension 
d = 3 of our universe is very special regarding the laws of 
gravity. This is the largest dimension of space at which 
all the series of equilibria of white dwarf stars (from clas- 
sical to ultra-relativistic) is stable. This may have impli- 
cations regarding the Anthropic Principle. 



III. SELF-GRAVITATING LANGEVIN 
PARTICLES 

A. The generalized Smoluchowski-Poisson system 

In this paper, we shall study a dynamical model of 
self-gravitating systems whose steady states reproduce 
the condition of hydrostatic equilibrium Eq. (w) . Specif- 
ically, wc consider the generalized Smoluchowski-Poisson 
system Q: 



dp 



(VP + pV$) 



A$ = S d G Pl 



(13) 



(14) 



where P(p) is a barotropic equation of state, i.e. the 
pressure P(r, t) depends only on the density of particles 
p(r, t) . This model describes a dissipative gas of self- 
gravitating Langevin particles in an overdamped limit 
£ — > +00 (where inertial effects are neglected) and in 



the thermodynamic limit N — > +00 (where the mean 
field approximation becomes exact) |l7|, Q. The GSP 
system is a particular example of generalized mean field 
Fokker-Planck equation [|5(J ■ It is associated to a stochas- 
tic process of the form 



dt e V pe 



R(i), 



(15) 



where R(i) is 
{Ri{t)R 3 (t')) = 



a white noise with (R(i)) = and 
SijS(t — if). This stochastic process de- 
scribes the evolution of each of the N Langevin particles 
interacting through the mean field potential $(r, t). For 
sake of generality, we have allowed the strength of the 
noise term in Eq. (15) to depend on the local distribu- 
tion of particles. This gives rise to anomalous diffusion 
and generalized pressure laws as discussed in |l9[ [50| . 

The Lyapunov functional (or generalized free energy) 
associated with the GSP system is 



F = J pf '^-dp'dr+^J p<S>dv. 



Easy calculations lead to 



F = - 



(VP + pV$) 2 dr < 0. 



(16) 



(17) 



The GSP system has the following properties: (i) the 
total mass is conserved, (ii) F < 0. (iii) F — <^> 
VP+pV^ = (hydrostatic equilibrium) <=> dtp = 0. (iv) 
p e? (r) is a steady state of the GSP system iff it is a criti- 
cal point of F[p] at fixed mass, (v) A steady state of the 
GSP system is linearly dynamically stable iff it is a (lo- 
cal) minimum of F[p] at fixed mass |74j. By Lyapunov's 
direct method |4£j , we know that if F~[p\ is bounded from 
below, the GSP system will relax towards a (local) min- 
imum of F[p] at fixed mass for t — > +00. If F[p] has 
several minima, the choice of the selected minimum will 
depend on a notion of basin of attraction: if the initial 
condition is sufficiently "close" to the minimum p eq (r), 
the distribution p(j,t) will converge towards p e q{^) for 
t —* +00. Finally, if F[p] has no global minimum (as can 
be the case for self-gravitating systems), the system can 
cither tend to a local minimum (mctastable) if it exists, 
or undergo collapse or evaporation. 

We are not claiming that this simple model accurately 
describes the dynamics of white dwarf stars or other as- 
trophysical systems. However, we have undertaken a sys- 
tematic study of the GSP system for different equations 
of state that have been considered in astrophysics. The 
main interest of this model is its simplicity (while being 
still very rich) which enables an accurate numerical and 
analytical treatment. This can be viewed as a first step 
before considering other, more realistic, dynamical mod- 
els of self-gravitating systems. On the other hand, in a 
completely different context, this model is isomorphic to 
the standard Keller-Segel model describ ing the chemo- 
taxis of bacterial populations (see Sec. VII). This is a 



further motivation to study this type of equations at a 
general level J75l . 



B. Isothermal spheres 



For an isothermal equation of state P = phsT/m, we 
recover the standard Smoluchowski-Poisson system pi[ : 



2r = v . 

dt 



l fh B T. 



-Vp + pV<& 



A$ = S d G P . 



(18) 



(19) 



Equation ( jl8| ) is an ordinary mean-field Fokker-Planck 
equation associated with a Langcvin dynamics of the 
form 



dr 
~dt 



-V$ 



l2k B T 
£m 



R(t), 



(20) 



where the strength of the noise is constant. The Lya- 
punov functional of the SP system can be written 



P 



F = k B T / — In — dr + - / p$ dr 



(21) 



This is the Boltzmann free energy Fb = E — TSb 
where E = (1/2) J p^ dr is the energy and Sb = 
—ks J(p/m)lii(p/m) dr is the Boltzmann entropy. The 
stationary solutions of the SP system are given by the 
Boltzmann distribution 



P 



AT*"* 



(22) 



where A is a constant determined by the mass M. These 
steady states can also be obtained by extrcmizing F at 
fixed mass, writing SF — a5M = 0, where a is a La- 
grange multiplier. The equilibrium distribution is ob- 
tained by substituting Eq. (22) into Eq. (19) leading to 
the Boltzmann-Poisson equation. Specializing on spher- 
ically symmetric distributions and defining 

p = Po e-^\ £ = r/r = (S d /3Gm Po )^ 2 r, (23) 

where po is t ne central density, we find after simple alge- 
bra that ip is solution of the Emden equation 



^- 1 d£ 



e 



d£ 



(24) 



with if) — and ip' — at £ = 0. The Emden equation can 
also be obtained from the fundamental equation of hy- 
drostatic equilibrium for an isothermal equation of state 
Ea, [Iq, pl| . Note that the isothermal spheres have a self- 
similar structure p(r)/po = e - ^ r / r °': if we rescale the 
central density and the radius appropriately, they have 
the same profile e~^^ . This property is called homology 

& 

For d = 1, the SP system is equivalent to the Burgers 
equation |23|, p0[ and it relaxes towards the Camm dis- 
tribution |6q | which is a global minimum of free energy 
for any temperature. For d > 2, there is no steady state 



with finite mass in an unbounded domain because the 
density of an isothermal self-gravitating system decreases 
as p ~ r~ 2 for r — * +oo [091. We shall thus enclose the 
system within a box of radius R |7q| . For box-confined 
systems, we must integrate the Emden equation (E4) un- 
til the normalized box radius £ = a with 



a = {S d pGmpo) 1/2 R. 



(25) 



It is useful to define a dimcnsionlcss control parameter 
j3GMm 



V = 



R d-2 



(26) 



Using the conservation of mass or the Gauss theorem, we 
get El: 



rj = aip'(a). 



(27) 



This equation relates the central density to the mass and 
the temperature. More precisely the relation rj(a) gives 
the mass M as a function of the central density (for a 
fixed temperature T) or the temperature T as a function 
of the density contrast K = p{0)/p(R) = e 4 '^ (for a 
fixed mass M). The curve 77(a) is plotted in Fig. 3 of |L2j. 
For 2 < d < 10, the series of equilibria 77(a) oscillates and 
presents a first turning point at r\ c = 77(01) (for d > 10, 
the series of equilibria does not display any oscillation) . 
According to Poincare's turning point argument |5j, p7| , 
configurations with a > a.\ are unstable (saddle points of 
free energy at fixed mass). This concerns in particular the 
singular isothermal sphere corresponding to a — > +00. 
Configurations with a < a\ are metastablc (local min- 
ima of free energy at fixed mass) and they exist only for 
V 5: Vc- There is no global minimum of free energy for 
self- gravitating isothermal spheres. For rj < rj c , depend- 
ing on the form of the initial density profile, the SP sys- 
tem can either relax towards a box-confined isothermal 
sphere (metastable) or collapse. This behavior has been 
illustrated numerically in Fig. 16 of |rj| . For -q > ?7 C 
the SP system undergoes gravitational collapse. This 
self-similar collapse, followed by thejormation of a Dirac 
peak, has been studied in detail in 
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If we remove 

the box, the SP system can either collapse or evaporate 
depending on the initial condi tion ( this behavior will be 
illustrated numerically in Sec. VIA ). 

The dimension d = 2 is critical and has been studied 
in detail in []12|, |l6| . The solution of the Emden equation 
is known analytically |68|]: 



,-■<!> - 



1 



1 



m 



(28) 



In an unbounded domain, the density profile extends to 
infinity but the total mass is finite because the density 
decreases as r -4 for r — ► +00. The total mass M = 



L p2nrdr is given by 



M = 



(3Gm 



,-i> 






where we have used the Emden equation (J2J) to get the 
last equality. Using Eq. (P8|), we find that £?/'' — ► 4 for 
£ — » +00. This yields a unique value of the mass (for a 
fixed temperature), or equivalently a unique value of the 
temperature (for a fixed mass) given by 



M r 



4k B T 
Gm ' 



knTr = 



GMr 



(30) 



For T = T c or M 
steady states 

p(r) 



M c , we have an infinite family of 



Po 



(l + |(r/r ) 2 



PotI 



k B T 
2nGm' 



(31) 



parameterized by the central density po- F° r Po ~~ * +00, 
we obtain a Dirac peak with mass M c . The steady states 
( pl| ) have the same value of the free energy, independently 
on the central density po (see Appendix 0) and they are 
marginally stable (8 2 F = 0). For T ^ T c or M ^ M c , 
there is no steady state in an infinite domain. For T > T c 
or M < M c , the solution of the SP system evaporates and 
for T < T c or M > M c , the solution of the SP system 
collapses. These different regimes have been discussed in 
detail in @,||. 

If we consider box confined configurations in d = 2, we 
observe that the control parameter (|2q ) is independent 
on the box radius and can be written 



M T 

n = BGMm = 4 — = 4— 

' M r T 



(32) 



Using Eqs. ( |27| ) and (p8[), we obtain the relation 77(a) = 
(a 2 /2)/(l + or/8) between the central density, the mass 
and the temperature. The density profiles are given by 
Eq. © with 8(r /i?) 2 = (T/T c - 1) = (M c /M - 1) 
so the central density is now determined by the mass M 
or the temperature T. Equilibrium states exist only for 

V < Vc = 4, i.e. M < M c or T > T c and, since the series 
of equilibria is monotonic, they are fully stable (global 
minima of free energy at fixed mass). In that case, the 
SP system tends to a box-confined isothermal sphere. For 

V = Vc = 4, i.e. M = M c or T = T C , the steady state is a 
Dirac peak containing all the mass. For 77 > ?7 C = 4 the 
SP system undergoes gravitational collapse (see Sec. 0). 

The mass-central density (for a fixed temperature) of 
two-dimensional isothermal spheres is plotted in Fig. g. 
We note the striking analogy with the mass-central den- 
sity of white dwarf stars in Fig. ||. Therefore, the critical 
mass ( p0| ) of isothermal spheres in two dimensions shares 
some resemblance with the Chandrasekhar mass. We 
shall show in the next section that this analogy (which 
is not obvious a priori) bears more significance than is 
apparent at first sight. 



C. Complete polytropes 

If we consider a polytropic equation of state P = Kp 1 
with 7 = 1 + 1/n, we get the polytropic Smoluchowski- 




FIG. 3: Mass as a function of the central density for two- 
dimensional box-confined self-gravitating isothermal spheres 
with fixed temperature. Equilibrium states exist only for 
M < M c . For M = M c , the density profile is a Dirac peak 
and for M > M c the system undergoes gravitational collapse. 
More precisely, this curve represents r/(a)/4 so it also gives the 
inverse temperature T c /T as a function of the density contrast 
1Z = po/p(R) = TZ(a) for a fixed mass. The corresponding 
density profiles are represented in Fig. 1 of pQ]. 



Poisson system pi[ : 
dp 



at 



= v 



{KVp 1 + pV$) 



A$ = S d Gp. 



(33) 



(34) 



Equation (|33|) is a generalized mean field Fokker-Planck 
equation associated with the stochastic process 



dr 

It 




2 -fp { ' 



where the strength of the noise depends on the local den- 
sity as a power law |]69| . The Lyapunov functional of the 
polytropic SP system can be written 



F 



K 



7-1 



(p>-p)dr + 



p$dr. 



(36) 



It can be interpreted as a generalized free energy of the 
form F = E - T eff S where E = (1/2) J>$dr is the 
energy, T e // = K is an effective temperature (polytropic 
temperature) and S = —1/(7 — l)/(p 7 — p) dr is the 
Tsallis entropy (the polytropic index 7 plays the role of 
the Tsallis q parameter). For 7 = 1, i.e. n — > +00, the 
polytropic equation of state P = Kp 1 reduces to P = 
Kp. It coincides with an isothermal equation of state 
P = pk B T jrn with temperature K = k B T ' jm leading to 
the standard Smoluchowski-Poisson system (p"8|)-(p"9|). 



The stationary solutions of the GSP system (t 
*iven by the Tsallis distributions 



7 



1 



1 V(T-l) 



K-f 



$ 



(37) 



where A is a constant determined by the mass M (by 
definition [x]+ = x if x > and [x]+ = if x < 0). These 
steady states can also be obtained by extremizing F at 
fixed mass, writing 5F—a5M = 0, where a is a Lagrange 
multiplier. The equilibrium distribution is obtained by 
substituting Eq. |}7|) into Eq. (||) leading to the Tsallis- 
Poisson equation. Specializing on spherically symmetric 
solutions and defining 



P = Pod n (0, Z = r/r 0l 



'u 



K{l + n) 

S d Gpl 1/n 



1/2 



(38) 

where po is the central density, we find after simple alge- 
bra that 9 is solution of the Lanc-Emden equation 






(39) 



with 6 = 1 and $' = at £ = 0. The Lane-Emden 
equation can equivalently be derived from the fundamen- 
tal equation of hydrostati c eq uilibrium with a polytropic 
equation of state |18|, |l9], 51 . Note that the polytropic 
spheres have a self-similar structure p{r)/po = 9 n (r/ro): 
if we rescale the central density and the radius appropri- 
ately, they have the same profile #"(£). This property is 
called homology Q ■ 

In this paper, we restrict ourselves to n > 0. Let us 
first discuss the case d > 2. For n > n^ = (d + 2)/(d — 2), 
unbounded self-gravitating polytropes have infinite mass 
because their density profile decreases like r~ a for r — > 
+00, with a = 2n/(n — 1). For n < rig = (d + 2)/{d — 
2), they are self-confined. In that case, the function 9 
vanishes at £ = £i and the density vanishes at _R„ = ro£i 
which defines the radius of the polytropc. The relation 
between the radius and the central density is 



i?* — 



A'(l + n) 
SdGp 



-i 1/2 



(40) 



r-R. 



The total mass M = L pSdr d x dr can be written as 



M = S dPo r% 



w s-d— 1 



e~ i dt - -s dP0 T%& 



dcd—lnl 



(41) 



where we have used the Lane-Emden equation ( |39| ) to 
get the last equality. Therefore, the relation between the 
mass and the central density is 



M = -S d po 



K{l + n) 

S d G P l- l/n 



d/2 



d-lnl 



£l 



(42) 



Eliminating the central density between Eqs. ( f40| ) and 
( [42] ) and introducing the index 



n 3 



d-2' 

we get the mass-radius relation 



(43) 



M (n-l)/n R [(d-2)(n 3 - n )]/n _ K(l + Tl) ^ (n _ 1)/n ,^. 



where 



GS 



_t(n+l)/(n-l)g/ 



l/i 



Let us introduce the polytropic temperature 



e 



K(l + n) 



iS 



l/n 



(45) 



(46) 



For < n < n 3 there is one, and only one, steady state 
for each mass M and temperature and it is fully stable 
(global minimum of F at fixed mass) . The GSP system 
will relax towards this complete polytrope (note that for 
n = 1 the radius i?„ of the polytropc is independent on 
the mass). For 713 < n < n§ there is one, and only one, 
steady state for each mass M and temperature O but it is 
unstable (saddle point of F at fixed mass) . In that case, 
the system will either collapse or evaporate. The index 
113 is critical. For n = 713, there exists steady solutions 
for a unique value of the mass (at fixed temperature 0): 



M r 



n 3 Q 
G 



"3/(13— 1) 



Ur, 



or for a unique temperature (at fixed mass M): 

(«3-l)/«3 



«3 V w "3 



(47) 



(48) 



M B 



(36/G) 



3/2 



U3 = 



For d = 3, we have 

10.487.. .{Q/Gf 2 and 6 C (G/3)(A//^ 3 ) 2/3 

0.20872. ..(G/M) 2 / 3 . As we have seen in Sec. || the 
Chandrasckhar limiting mass of relativistic white dwarf 
stars is connected to the limiting mass (|47J ) of critical 
polytropes. For a polytropic equation of state with crit- 
ical index n = 713, and for M = M c , we get an infinite 
family of steady solutions 



p(r)=p e n *(r/r ), 



Po4 



1 
~Sd 



en 3 y/ 2 



G 



(49) 



parameterized by the central density po- For po —* +00, 
the density profile tends to a Dirac peak with mass M c . 
These solutions have the same equilibrium free energy 
F[Peq] = —dKM/(d — 2) independently on the central 
density po (see Appendix [A]) and they are marginally 
stable (S 2 F = 0). For M < M c (at fixed temperature) or 
> 6 C (at fixed mass), the solutions of the GSP system 



evaporate and for M > M c (at fixed temperature) or 
8 < O c (at fixed mass), they collapse. These different 
regimes will be studied in detail in Sees. |v| and VI. 

For d = 2, we find that 713 — > +00, so we realize that 
isothermal systems (n = +00) in two dimensions are sim- 
ilar to critical polytropes (n = 77,3) in higher dimensions 
d > 2. This is why the critical mass of isothermal spheres 
in d = 2 shares some analogies with the Chandrasekhar 
mass in d = 3 since they both correspond to critical poly- 
tropes with index n — n 3 plj ]. Comparing Eq. (29) with 
Eq. ( pl|) we find that for d — > 2 and n — n 3 — > +00, we 
have the limit 



lim ti 3 w„ 3 



4. 



(50) 



This limit can also be obtained from Eq. (79) of |pl| . 
With this relation, we find that the critical mass and 
the critical temperature in d = 2 given by Eq. (J3CJ) are 
particular cases of Eqs. (R^j) and (Eq). 

Finally, for d = 1 with n > (and for d = 2 with 
< n < +00), the GSP system always relaxes towards 
a complete polytrope which is a global minimum of free 
energy. Thus there is no critical dynamics for d < 2 (and 
for d — 2 with n 7^ +00). 



is valid only for incomplete polytropes whose density pro- 
file is arrested by the box (i.e. p(R) > 0). For n > ns, 
this is always the case. For < n < 775 , using the identity 



£1 



R 

R* 



(55) 



the polytrope is confined by the box if i?„ > R, i.e. a < 
£1. For i?» < i?, i.e. a > £1, we have complete polytropes 
whose density profile vanishes before the wall. In that 
case, we need to integrate the Lane-Emden equation until 
the natural polytropic radius £ = £1. For a > £1, the 
relation (K4) is replaced by 



'/ 



,n/(n-l). 
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(d-2)(n-n 3 )/(n-l) 



(a > £l), (56) 



which is equivalent to the mass-radius relation (Q). Us- 
ing Eq. (p5J) , it can be expressed in terms of a, giving 
the relation between the mass and the central density (at 
fixed temperature) for complete polytropes. Finally, the 
intermediate case is i?* = R, i.e. a = £1, at which the 
density profile vanishes precisely at the box radius. In 
that case, we have 



>l 



n/("-l), 



(" = 6)- 



(57) 



D. Box confined polytropes 

For systems confined within a box of radius R, we need 
to integrate the Lane-Emden equation (|3^) until the nor- 
malized box radius £ = a with 



Q' = R/tq = 



SdGp 



1-1/fl 



isT(n+l) 



1/2 



i?. 



(51) 



It is useful to define a dimcnsionless control parame- 
ter (the definition of this parameter has been slightly 
changed with respect to our previous paper [pi]): 



7? = 71/ 



#(l + n) 



-1 n/(n-l) 



1 



R (d-2)(n-n 3 )/(n-l) 



(52) 



In terms of the polytropic temperature (|46|), it can be 
rewritten 



V 



Qn,/{n-X) M 



Qn/(n-l) R (d-2)(n-n a )/(n-l) ' 



(53) 



Note that for n — > +00, we have = if = /cbT/to 
and the definitions (pO) and (53) coincide. Using the 
conservation of mass or the Gauss theorem, we get [plj : 



V 



= -n'^^-^a^+^/^-^B'io), (a < &). (54) 



This equation relates the central density to the mass (at 
fixed temperature and box radius). In fact, this relation 
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FIG. 4: Series of equilibria for box-confined polytropes with 
different index (the figure is done for d — 3). The full lines 
(a. < £1) correspond to incomplete polytropes whose profile is 
arrested by the box and the dashed lines (a > £1) correspond 
to complete polytropes that are self-confined. 

The relation 77(a) defines the series of equilibria con- 
taining incomplete (for a < £1) and complete (for a > £1) 
polytropes. It gives the mass M as a function of the cen- 
tral density (for a fixed temperature and box radius R) 
or the temperature as a function of the density contrast 
(for a fixed mass M and radius R) . Different examples of 
curves 77(a) are represented in Fig. |4| for various indices 
in d = 3: 

• For n < 77,3, the series of equilibria 77(a) is monotonic. 
Since polytropic spheres are stable in absence of gravity 



10 



(corresponding to a — ► 0) and since there is no turning 
point, the Poincare argument implies that all the poly- 
tropes are stable. It can be shown furthermore that they 
are fully stable (global minima of free energy at fixed 
mass) so that the GSP system will tend to a steady state 
for t — > +00. For r\ < 771 = 77(^1) = ?W("~ 1 )cj„, the 
GSP system tends to an incomplete polytrope confined 
by the box. For r\ > 771, the GSP tends to a complete 
polytrope with radius i?„ < R. This has been illustrated 
numerically in Fig. 21 of Q for n = 3/2 in d = 3. This 
index corresponds to a classical white dwarf star in as- 
trophysics. If we remove the box, the GSP system always 
tends to the complete polytrope. 

• For n > 713, the series of equilibria 77(a) presents a 
turning point at rj c = 77(01). According to the Poincare 
turning point argument, configurations with a > u\ are 
unstable (saddle points of free energy at fixed mass). 
This concerns in particular the case of complete poly- 
tropes for 713 < 77 < 77,5 (corresponding to a = £1), the 
Schuster polytrope n = 77,5 and the singular polytropic 
spheres for n > n 5 (corresponding to a — +00). Con- 
figurations with a < a\ are metastable (local minima of 
free energy at fixed mass) and they exist only for 77 < r/ c . 
There is no global minimum of free energy for 77 > 713. 
For 77 < 77 c , depending on the form of the initial density 
profile, the GSP system can either relax towards an in- 
complete polytrope confined by the box (metastable) or 
collapse. For 77 > ?7 C , the GSP system undergoes gravita- 
tional collapse. This self-similar collapse has been stud- 
ied in detail in |pl|| . It is very similar to the self-similar 
collapse of isothermal systems in d > 2 corresponding to 
77 — > +00. If we remove the box, the GSP system can 
either collapse or evaporate depending on the initial con- 
dition (this will be illustrated numerically in Sec. [Vj). 

• The case 77 = 773 is critical and will be studied in 
detail in this paper. For the critical index 77 = 773, the 
control parameter is independent on the box radius and 
can be written 



'/ 



M 



ri3/("3 — 1) 



(58) 



In terms of the critical mass 
(Eq), we have 



]) or critical temperature 



n 3 /(n3-l) 



M 



77" 3/{ " 3 ~V 



"3/(13-1) 



(59) 

For incomplete polytropes with a < £1, the relation 77(a) 
between the central density, the mass and the tempera- 
ture is given by Eq. (p4h. Their density profile is given 
by Eq. ( |49| ) where r is determined by (0 c /0) rf / 2 = 
M/M c = -(l/u} n3 )(R/r ) d ~ 1 d'{R/r ), equivalent to re- 
lation (|5J), so the central density is now determined by 
the mass M or the temperature 0. Complete polytropes 
with a > £1 exist for a unique value of the control pa- 
rameter 



This corresponds to the critical mass M = M c or critical 
temperature = C . Equilibrium states exist only for 
V < rjc, i-c M < M c or > C . For 77 < ?7 C , they 
are fully stable (global minima of free energy at fixed 
mass). In that case, the GSP system relaxes towards an 
incomplete polytrope confined by the box. For 77 = 77c, i.e 
M = M c or = C , we have an infinite family of steady 
states parameterized by their central density a > £1 or 
cquivalcntly by their radius i?* < R. They arc marginally 
stable {S 2 F = 0). For 77 > ?7 C , i.e M > M c or < C , 
the GSP system undergoes gravitational collapse. The 
collapse dynamics is expected to be similar to the critical 
collapse of isothermal systems with 77 — > +00 in d = 2 
(see below). If we remove the box, the solution of the 
GSP system evaporates for r/ < rj c , i.e. M < M c or 
> C and collapses for 77 > 77,-, i.e. for M > M c or 
< C . These different regimes will be studied in detail 
in Sees. [v]and[vi. 
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FIG. 5: Mass as a function of the central density for box- 
confined self-gravitating polytropic spheres with critical index 
71 = 773 = 3 in d — 3. Incomplete polytropes with p(R) > are 
represented by a solid line and complete polytropes with R, < 
R are represented by a dashed line. For po —> +00, the density 
profile tends to a Dirac peak. Equilibrium states exist only for 
M < M c - For M > M c the system undergoes gravitational 
collapse. The curve represents 77(a)/[7i3 33 ~ lj„ 3 ] so it also 
gives the inverse temperature (0 C /O)™ 3 ' '™ 3-1 ' as a function 
of the density contrast lZ(a) for a fixed mass. 

The mass-central density relation (for a fixed temper- 
ature) of box-confined self-gravitating polytropic spheres 
with critical index n = 77,3 is plotted in Fig. |5j and the 
corresponding density profiles (illustrating the notion of 
complete and incomplete polytropes) arc plotted in Fig. 
0. We note the striking analogy with the mass-central 
density relation of white dwarf stars in Fig. [2| Indeed, 
ultra-relativistic white dwarf stars are equivalent to poly- 
tropes with critical index 77, = 77,3 = 3 in d = 3. In this 
context, the critical mass M c corresponds to the Chan- 
drasekhar limit. We emphasize, however, that we are 
considering here pure critical polytropes enclosed within 
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FIG. 6: Density profiles of complete and incomplete poly- 
tropes for the critical index n3 = 3 in d = 3. We have consid- 
ered three values of the central density po = (M c /SdR d Un 3 )a d 
corresponding to a — 3 < £i (incomplete polytrope: i?» > R, 
M < M c ), a = £l = 6.89685... (limit polytrope: i?» = i?, 
M = M c ), and a = 20 > £i (complete polytrope: R, < i?, 
M = M c ). 




FIG. 7: Mass as a function of the central density for crit- 
ical polytropes n — U3 for different dimensions of space. 
We have plotted r/ — n 3 uj n3 M/M c as a function of a = 
(poSdR o)n 3 /M c ) • The maximum mass is reached at the 
bullet corresponding to a = £1 , 77 = r\ c — n 3 u„ 3 . 



a box while in Sec. || we considered self-confined par- 
tially relativistic white dwarf stars for which a box is 
not needed. It is only when M — > Mchandra (ultra- 
relativistic limit) that they become equivalent to pure 
polytropes. Furthermore, at M = Mchandra for white 
dwarf stars, the only steady state is a Dirac peak while 
at M = M c for pure critical polytropes, we have an infi- 
nite family of steady states with different central densi- 
ties (the same difference holds between critical polytropes 
n = i%3 in d > 2 and isothermal spheres n — n^ = +00 
in d = 2; compare Figs, ra and ph. Finally, in Fig. M, 
we plot the mass as a function of the central density for 
different dimensions of space d. This figure illustrates 
in particular the connexion between the critical mass in 
d = 3 reached for a finite value of the central density and 
the critical mass in d = 2 reached for an infinite value of 
the central density. 



IV. THE CRITICAL INDEX FROM 
DIMENSIONAL ANALYSIS 

It is instructive to understand the origin of the critical 
index 74/3 = 2(d— l)/2 or 77.3 = d/(d— 2) from simple di- 
mensional analysis. Here we consider unconfined systems 
in d dimensions with arbitrary value of 7. The polytropic 
Smoluchowski-Poisson system can be written 



dp 
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K'ypi-^Vp + pV$) 



A$ = S d G P . 



-V-J, (61) 



(62) 



3d = —K'fp' 1 ~ 1 yp and a gravitational drift J 9 = — pV$. 
Based on dimensional analysis, the diffusion current can 
be estimated by 

J d ~ +K 1 (M/L d )->- 1 (p/L) ~ +(l/i)^- 1 )+ 1 , (63) 



and the drift term by 

J g = -pGM/L*- 1 ~ -(l/L) d -\ 



(64) 



where M is the mass of the system and L is the charac- 
teristic size of the system. 

The system will collapse to a point if gravity overcomes 
(anomalous) diffusion, i.e. \J g \ ^S> \Jd\-, when L — ► 0. 
This will be the case if d— 1 > ^(7— 1) + 1, i.e. 7 < 74/3. 
Conversely, if 7 > 74/3, the diffusion term can stabi- 
lize the system against gravitational collapse so that the 
system can be in stable equilibrium. The system will 
evaporate to infinity if (anomalous) diffusion overcomes 
gravity, i.e. \J&\ 3> \J g \, when L — > +00. This will be the 
case if ^(7— 1) + 1 < d— 1, i.e. if 7 < 74/3. Conversely, if 
7 > 74/3; the gravitational attraction can prevent evap- 
oration so that the system can be in stable equilibrium. 
In conclusion, we find that the system can be in a stable 
equilibrium state iff 7 > 74/3, i.e. 1/rc > l/n^. In the 
opposite case, the system can either collapse to a point 
or evaporate to infinity. By this very simple argument, 
we recover the stability criterion of self-gravitating poly- 
tropic spheres obtained by other methods (see Appendix 
B of 0). 

The critical case is obtained when Jd ~ J g implying 
d("f— 1) + 1 = d— 1, i.e. 7 = 74/3 or, cquivalently, n = n^. 
In that case, the stability of the system will depend on 
its mass. The system will collapse to a point if gravity 



The current J = J d + J 3 appearing in the Smoluchowski overcomes diffusion, i.e. \J g \ > \J d \, when L — > 0. This 
equation is the sum of two terms: a diffusion current will be the case if M > M c , where M c ~ (K/G) d / 2 is a 
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critical mass. The system will evaporate to infinity (in an 
unbounded domain) if (anomalous) diffusion overcomes 
gravity, i.e. \J d \ » \J g \, when L — ► +00. This will be 
the case if M < M c . Therefore, at the critical index 
7 = 74/3 i- e - n = ns, the system collapses if M > M c 
and evaporates if M < M c . Again, this is fully consistent 
with the results obtained in Appendix B of pW. 



V. COLLAPSE DYNAMICS 

For < n < 713 in a space with d > 2 dimensions, the 
GSP system tends to an equilibrium state. For n > n 3 , 
it can undergo gravitational collapse. For n > n 3 with 
d > 2, the collapse is self-similar as studied in |5l| (the 
case of negative indices n < is studied in [p2f j. In 
the present section, we consider the collapse dynamics 
of self-gravitating Langevin particles associated with the 
critical index 713 = d/(d — 2) in d > 2 dimensions which 
presents non trivial features. 



It is also convenient to introduce the function s(r,t) = 
M(r,t)/r d which has the same dimension as the density 
and which satisfies 

ds_ G f ds \ 1/n fd 2 s d+l ds\ f ds \ 
dt \ dr J \dr 2 r dr J \ dr J 

(70) 
For n — > +00, these equations reduce to those studied in 
Refs. [[n], [12| in the isothermal case. 

When the system collapses, it is natural to look for 
self-similar solutions of the form 



p(r,t) = Po (t)f 



r»(t)J' 



ro = 



e 



1-1/71 

Po 



1/2 



■ (71) 



The relation between the core radius ro and po (propor- 
tional to the central density ||77| ) is obtained by requiring 
that the diffusive term and the drift term in Eq. ( pq ) 
scale in the same way. This relation can be rewritten 
Po r o ~ 1 with 



A. Generalities: self-similar analysis 

From now on, we adopt normalized variables such that 
G = M = R = £=1. The unique control parameter is 
the temperature 8. For spherically symmetric solutions, 
using the Gauss theorem, the GSP system can be written 
in the form of an integrodiffcrcntial equation 



dp 



1 d 

r d—l Q r 



„d-l 



{S d pf' n @^ 



r d-l 



p(r')S d r ld - 1 dr 



Introducing the mass within a sphere of radius r 



M(r,t)= / P {r')S d r' d - 1 dr' , 
Jo 



(65) 



(66) 



the GSP system can be formulated through a unique non- 
linear dynamical equation for M(r,t): 



dM 



e 



1 dM 

r d-l Q r 

M dM 



r d 1 Q r 



!/« f^2 



d 2 M d-ldM 



dr 2 



dr 



(67) 



If the system of total mass M = 1 is confined within a box 
of radius R — 1, the appropriate boundary conditions are 



M(0,t) = 0, M(l,i) = l. 



(68) 



If the system is not confined, the second condition should 
be replaced by 



M(oo,£) = 1. 



(69) 



2n 

a = . 

n- 1 

In terms of the mass profile, we have 



(72) 



M(r, t) = M (t)g (-^ j , with M Q (t) = p Q r d , 

(73) 



and 



g(x)= / f{x')S d x' d - x dx'. 
Jo 

In terms of the function s, we have 



(74) 



with S(x) = 9 M 
r (t) J x a 



s(r,t) = p (t)S 



Inserting the ansatz ( |75| ) in Eq. ( \I(\ ) and using Eq. J7J] ) , 
we obtain 



15) 



1 dpp _ 
p 2 dt 



(76) 



and 



aS + xS' = {xS' + dS) 1/n (s" + ^-S') + (xS' + dS)S. 

_ (77) 

Assuming that Eq. ( |77|) has a solution so that the self- 
similar solution exists, Eq. (W&j is readily integrated in 



p (t) = -{tcoll -t) \ 



(78) 



implying a finite time singularity. On the other hand, the 
invariant profile has the asymptotic behavior f(x) ~ x~~ a 
for x — > +00. 
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B. The two-dimensional isothermal case 

In d = 2 dimensions, the critical index is 713 = +00 
corresponding to the isothermal case studied in [12 ] (in 
that case 9 = T). Since the study of the critical dynam- 
ics is rather complicated, it can be useful to summarize 
our results, with some complements and amplifications, 
before treating the case d > 2. 

In d = 2, there exists a critical temperature T c = 1/4. 
If the system is enclosed within a box and T > T c , it re- 
laxes to an equilibrium distribution confined by the box. 
If the system is not confined and T > T c , an evapo- 
ration process develops which has been studied in [p.6| . 
For T = T c , the system undergoes gravitational collapse. 
The evolution is self-similar and leads to a Dirac peak 
containing the whole mass M = 1 for t — ► +00. In a 
bounded domain, the central density grows exponentially 
p2| rapidly with time and in an unbounded domain, the 
central density increases logarithmically pa with time 
(and a tiny fraction of mass is ejected at large distances 
to satisfy the moment of inertia constraint at T = T c ). 
Note that the Dirac peak is also the stationary solution 
of the SP system at T = T c . 

For T < T c , and irrespectively of the presence of a 
confining box, there is no steady state and the system 
collapses. Looking for an exact self-similar solution of 
the form ( |7l| ) we obtain par 2 — T, a — 2 = d and a 
scaling equation 



S" + -S') + (xS' + 2S)(S 
x 



1 



0. 



(79) 



However, this equation does not have any physical solu- 
tion for large x. In fact, this could have been anticipated 
from the fact that the scaling functions s(x) and f{x) 
should decay as x~ 2 = x~ d for large x. Then, the to- 
tal mass in the profile is of order po r o / ° x~ 2 xdx ~ 
ln(l/ro), which unphysically diverges when ro goes to 
zero. Said differently, the scaling profile at t = t co ii is 
p ex r~ 2 so that the mass M = J p(r)2irrdr diverges log- 
arithmically for r — > 0. This logarithmic divergence is 
symptomatic of the formation of a Dirac peak resulting 
from a pseudo self-similar collapse. In the case d = 2, this 
situation can be analyzed analytically in great detail. 

To that purpose, we note that the profile which cancels 
out the r.h.s. of the SP system is exactly given by 



M 1 (r,t)=AT 



(r/r (t)) 2 
l + (r/r (t)y 



Pi{r,t) 



4po(*) 



1 



it (l + (r/r (i)) 2 ) 2 



with 



Po{t)r {t) 2 =T. 



(80) 



(81) 



(82) 



If we consider time independent solutions (dp/dt = 0) 
and impose the conservation of mass, we recover the 



steady solutions which exist for T > T c in a bounded 
domain (in that case r = {T/T c - l) 1 ' 2 ) and for T = T C 
only in an infinite domain (in that case we get a family 
of distributions parameterized by ro). However, in the 
present case, we consider the case T < T c and seek the 
temporal evolution of po(t) and ro(t). We argue that the 
solution (pll) gives the leading contribution of the den- 
sity profile in the core. This profile contains a mass T/T c . 
We expect that the collapse will lead to po(t) — ► +00 and 
ro(t) — > for t — > t co u (finite time singularity). Then, we 
see that the profile (^l|) leads to a Dirac peak with mass 
T/T c , i.e. 



Pi(M) 



T 
T, 



5{v). 



(83) 



The excess of mass will be contained in the profile ex- 
tending in the halo. Therefore, we look for solutions of 
the form 

p(r,t) = pi(r,t) + p 2 (r,t), 

= Po(t)fi(r/r (t)) + p (tr^/ 2 h(r/r (t)). 

(84) 

The first component has a scaling behavior and domi- 
nates in the center of the collapse region. It leads to a 
Dirac peak containing a fraction M c = T/T c of the total 
mass M = 1 at t = t C oll- The second component obeys 
a pseudo-scaling and J2{x) ~ x~ a ^ for large x, with 
an effective scaling exponent a(t) which very slowly ap- 
proaches the value 2 (expected from the naive self-similar 
analysis) when t — > t co u. Thus, at t = t co u, we get 



p(r,t) ^ M c S{r) + X (r,t), 



(85) 



where yfr) is singular at r = behaving roughly as r~ 2 . 
In Fig. g, we illustrate this decomposition of the density 
profile into two components. It is shown in |L2) that the 
central density satisfies an equation of the form 



1 dpo att)/2 

Po dt 



(86) 



instead of Eq. J7q), and that the effective scaling expo- 
nent a(t) depends on the central density as 



e(t) = 1 



a{i) /lnlnpo(^) 



2 y 2\np {t) 

This yields po ~ (t co u — t)~ 1+c ^ or equivalently 



(87) 



„, / / rs l\ In r I In llnrl 
ln(p r) ~ -21n(r /V^) ~ ^ l - L_J L, 



(88) 



where we have noted r = t co u — t. 

Prior to our work |12J , and unknown to us at that 
time, Herrero & Velazquez pOf had investigated the same 
problem in the context of chemotaxis using a different 
method based on match asymptotics. For T < T c (as 
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C. The critical polytropic case with d > 2 



10 



lO" 2 io 1 

r 



10 



We now consider the critical index n = iij, = d/(d— 2) 
with d > 2. There exists a critical temperature C = 
l/[n 3 J r ^ 1)/n3 ] (in d = 3, we have C = 0.20872...). 
If the system is confined within a box and O > O c , it 
relaxes to an incomplete polytropc. This is illustrated 
in Fig. |9[ If the system is not confined and 9 > C , 
an evaporation process develops which will be studied in 
the next section. In the confined case, when the general- 
ized temperature reaches the value C , the equilibrium 
density profile vanishes exactly at R = 1. For < C , 
and irrespectively of the presence of a confining box, the 
system collapses. 



FIG. 8: For d — 2, n — ns =+oo, and deep into the collapse 
regime for T = T c /2 = 1/8, we plot the density profile (full 
line), emphasizing its two components: the core is dominated 
by the invariant scaling profile (dotted line) given analytically 
by Eq. (plh containing a mass M c = T/T c , and the halo obeys 
pseudo-scaling (dashed line) with an exponent a(t) tending 
slowly to d = 2 as po — » +co. 



far as we know, they did not consider the case T = T c 
treated in |l2|), they showed that the system forms a 
Dirac peak of mass M c = T/T c (within our notations) 
surrounded by a halo containing the excess of mass. On a 
qualitative point of view, the two scenarii are consistent. 
On a quantitative point of view, however, the scaling laws 




ln(p r) ~ -21n(r /v^) ~ ^2| ln-r| 



V\^' 



In I lnrl 



(89) 



obtained by Herrero & Velazquez (HV) are slightly dif- 
ferent from ours (SC). They lead to an effective exponent 
given by 



1- 



a(t) 



lnpo 



Ul--±=)^, (90) 
2 V Vlnpo/ ln Po 



instead of Eq. (|87|). For the densities accessible nu- 
merically, one gets ascipo = 10 3 ) = 1.252... while 
a H v(po = 10 3 ) = 0.751... and a sc (po = 10 5 ) = 1.348... 
while cxhv(Po = 10 5 ) = 1.017.... Numerical simulations 
performed in p2| show a good agreement with the pre- 
dicted values of asc for the densities accessible. How- 
ever, in view of the complexity of the problem, and of the 
logarithmic (and sub- logarithmic!) corrections, it is diffi- 
cult to understand the origin of the (slight) discrepancy 
between the two approaches. In any case, they both show 
that the collapse is not exactly self-similar but that the 
apparent scaling exponent a(t) is a very slowly varying 
function of the central density. 



FIG. 9: In d = 3 and for n 3 = 3 and 6 = 0.25 > 9 C in a 
finite box (R — 1), we show the density at successive times, 
illustrating the convergence to the equilibrium density profile 
(dashed line). The insert illustrates the exponentially fast 
saturation of the central density for O > C , whereas a slower 
algebraic saturation is expected right at Q = G c . 

We can naively look for self-similar solutions of the 
form described in Sec. 



VA. For n = n3, we find a = d, 



Po r o 



Q d / 2 and the scaling equation 
S" + *±ls> + (xS' + dSf/ d (S - 1) = 0. 



(91) 



It happens that as in the case (d = 2, 713 = 00), this 
equation does not have any physical solution for large 
x. Again, this could have been anticipated from the fact 
that the scaling functions s{x) and f(x) should decay as 
x -^n 3 /(n 3 -i) _ x -d^ f or i ar g C x _ Then, the total mass in 

the profile is of order 



Po^o 



l/ro 



- d xx d - 



1 dx ~ ln(l/r ) 



(92) 



which unphysically diverges when r$ goes to zero. Said 
differently, the scaling profile at t = t co u is p oc r~ d so 
that the mass M = J p(r)Sdr d ~ 1 dr diverges logarithmi- 
cally for r -> f7l. 
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FIG. 10: For d = 3, n = nz = 3, and deep into the collapse 
regime for O = 0.75S C , we plot the density profile (full line), 
emphasizing its two components: the core is dominated by the 
bounded invariant scaling profile (complete polytrope of index 
723) containing a mass M c — (0/O c ) (dotted line), and the 
halo obeys pseudo-scaling (dashed line) with an exponent a(t) 
tending slowly to d = 3 as po — » +00. 



Hence, for n = 713 in d > 2, we expect a situation sim- 
ilar to what was obtained for (d = 2, 113 = 00). However, 
the situation is more difficult to analyze because the sta- 
tionary profile is not known analytically in the present 
case (this analytical profile was at the basis of our anal- 
ysis in [|12|). Using the results of Sec. Ill C , the profile 
which cancels out the r.h.s. of the GSP system is given 
by 



d/2 



Pi(r,t) 



S,, 



■p (t)e^(r/r (t)), 



with 



Po(t)r (t) d = Q d ' 2 



(93) 



(94) 



If we consider time independent solutions (dp/dt = 
0) and impose the conservation of mass, we recover 
the steady solutions that exist for > 8 C in a 
bounded domain (in that case, we have (9 c /9) d / 2 = 
-(l/u n3 )(R/r ) d - 1 6'(R/r )) and for 9 = 9 C only in an 
infinite domain (in that case we get a family of distri- 
butions parameterized by ro). However, in the present 
case, we consider the case 9 < 9 C and seek the temporal 
evolution of po(t) and ro(t). We argue that the solution 
( P3| ) gives the leading contribution of the density profile 
in the core. This profile vanishes at R*(r) = £iro(£), has 



a central density (n 3 
Sec. pICl): 



a/ 2 



/S c i)pa(t) and contains a mass (see 



Mr 



e \ d/2 

97 



(95) 



We expect that collapse will lead to p a {t) — > +00 and 
ro(t) — > 0. Then, we see that the profile (|93[) tends to a 



Dirac peak with mass M c , i.e. 



(-) 



MM)-(q- 



c//2 



5{v) 



(96) 



The excess of mass will be contained in the profile ex- 
tending in the halo. Therefore, we look for solutions of 
the form 

p(r,t) = pi{r,t) +p 2 (r,t), 

= Po(t)h(r/r (t))+p (tr^/ d Mr/r (t)). 

(97) 

The first component [|79| has a scaling behavior and dom- 
inates in the center of the collapse region. It leads to a 
Dirac peak containing a fraction M c = (9/9 c ) d / 2 of the 
total mass at t — t co u. The second component obeys 
a pseudo-scaling and f<z(x) ~ x~ a ^ for large x, with 
an effective scaling exponent a(t) which very slowly ap- 
proaches the value d (expected from the naive self-similar 
analysis) when t — > t co u . At t — t co ii , the first component 
p\ (r, t) tends to a Dirac peak at the origin containing the 
mass M c , whereas the second component develops a sin- 
gularity at r — 0. Thus, we have 

p(r,t) ^ M c d(v) + X (r,t), (98) 

with x( r ) behaving roughly as r~ d . In Fig. |l0|, we illus- 
trate this decomposition of the density profile into two 
components. 
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FIG. 11: For 6 = 0.756 c (here d = 3, n = n 3 = 3), we 
plot Pq 1 ^)^- (top full line) and po(t) (bottom full line) as 
a function of po(t)- Both grow with an effective exponent 
a/3 ~ 0.93 (dotted lines), which slowly increases and should 
saturate to unity (the dashed line has slope unity). 

In Fig. [Ill, we show that perfect scaling which would 
imply PQ l {t)-^- ~ Po is not obeyed. Instead, in the 
accessible density range, PQ l {i)-^- decays with an ap- 
parent power-law of po which increases very slowly with 
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time, but remains less than unity. We expect to have a 
relation of the form 



dimension of 0: 



1 dp a(t)/d 

Po at 



(99) 



which is indeed confirmed by the numerics. In Fig. O, we 
also plot the central density in the pseudo-scaling com- 
ponent 



Po(t) = P a a (t)/d (t), 



(100) 



which shows that the effective exponent a(t) slowly con- 
verges to a = d. 




FIG. 12: For 9 = 0.759 c (here d = 3, n = n 3 = 3), we plot 
/j ( x ) ( as defined in the text) as a function of x = r/ro, for 
different times for which the central density evolves from 10 
to 10 7 . Pseudo-scaling is observed. The envelop of the tails 
decays with an apparent exponent a w 2.8 (right dashed line), 
while the small x behavior is quadratic (left dashed line). 

Finally in Fig. [12], we display the apparent scaling be- 
havior of p2(r, t) = po(t) a ^ t ^ d f2(r/ro(t)), associated to a 
value ofaK 2.8, fully compatible with the value obtained 
in Fig. ^ (in d = 3). 



VI. EVAPORATION DYNAMICS IN 
UNBOUNDED SPACE 

A. The case n > 713 

When the system is not confined to a finite box, the 
nature of the dynamics crucially depends on the value 
of the polytropic index n with respect to n^. As before, 
we consider d > 2 and n . > 0. If n < 77,3, there exists 
equilibrium solutions (fully stable complete polytropes) 
which are reached for any initial density profile. If n > 
n 3 , depending on the initial density profile and on the 
temperature, the system can collapse or evaporate. If 
Ro is the typical extension of the initial density profile 
containing a mass M, one can form a quantity with the 



e* 



Q M {n-V>/n 



R 



(d-2)(n-n 3 )/n ; 



(101) 



which plays the role of an effective critical temperature. 
If <C 0*, the system should collapse as it would do if 
confined in a box of typical radius Ro [pl| . If 3> 0*, 
the system should evaporate in the absence of an actual 
confining box. Hence, for a given initial profile, there 
exists a non universal 0* separating these two regimes. 
We present numerical simulations for the case n > 773 . In 
Fig. [13, and for a particular initial process, we illustrate 
the fact that depending on the value of with respect to 
a non universal 0*, the system can collapse or evaporate. 
In the evaporation regime and for n > n 3 , a scaling anal- 
ysis shows that gravity becomes gradually irrelevant and 
that this process becomes exclusively controlled by free 
(anomalous) diffusion. This fact is illustrated in Fig. [lj. 
Indeed, when the evaporation length ro(t) — ► +00, we 
see from Eq. (g7() that the gravitational term becomes 
negligible in front of the diffusion term: 



M dM 

r d—\ Q r 



<0 



1 dM 



Or 



-d-l 



1/n 



8 2 M 

dr 2 



(102) 



if d > d/n + 2, i.e. n > n 3 . Therefore, for t > 1, the GSP 
system reduces to the pure anomalous diffusion equation 



dp 

Of 



K d 

-*(Z — X {"if* 



,d-X d _f_ 

dr 



(103) 



with K = Sj~ 0/7. This equation has self-similar so- 
lutions that were first discovered by Barenblatt [|70[ in 
the context of porous media. These solutions are closely 
related to the form of generalized thermodynamics intro- 
duced by Tsallis || . 

Using the original idea of Pl astino & Plastino WW , we 
look for solution of Eq. (102) in the form of a Tsallis 
distribution with index 7 and time dependent coefficients 



1 



p(r,t) = -p (t)[l-( 1 



l)(r/ro(*)) 2 ]i /(7 " 



(104) 



For 7 > 1, i.e. n > 0, we have a profile with com- 
pact support where the density vanishes at r max (t) = 
r o(t)/Vl ~ 1- F° r 7 < 1, i.e. n < 0, the density decreases 
like p ~ r" 2 ^ 1 r ' and the total mass is finite provided 
that 7 > 71/3 = (d — 2)/d, i.e. n < —d/2. Requiring 
that the profile (104) contains all the mass M = 1, and 
imposing 



Po{t)r (t) d 



1, 



(105) 



with find the normalization factor 




Z = 



[1 - (7 - l^^^SdX^dx. (106) 



Then , su bstituting the ansatz (104) with Eq. ( |105[ ) in 
Eq. (103), we obtain 



Po 



-2ds 1 d - l ez l -'<pl +2,d . 



(107) 
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FIG. 13: In d — 3 and n — 5 > 713, and for a given initial den- 
sity profile (M(j-) = r /(e _r + r 2 ) ' ; fat line), we show the 
collapse dynamics observed at = 0.15 (full lines for different 
times before t co ll) and the evaporation dynamics observed at 
= 1 (dashed lines for different times). For this particular 
initial condition, we find O* ~ 0.206. 




FIG. 14: In d = 3 and n = 5 > 713, we present the evaporation 
density data collapse at O = 1. As time proceeds, the effect of 
gravity becomes negligible, and the scaling profiles converge 
to the one corresponding to free diffusive evaporation (full 
line). This is the Barenblatt so lutio n whose invariant profile 
is a Tsallis distribution of Eq. (114) with index 7. 



Solving this equation with the initial condition p(r, t = 
0) = <5(r), we get 

1 



Po(t) 



[2d(7 - ii /3 )sj- ez'-n] 1 /^-^) 



(108) 



This is valid for 7 > 71/3, i.e. n > or n < —d/2. We 
note the scaling laws for large times: 



Po 



(t) ~ r dn/(d+2 " ) , ro(t) ~ t n l (d+2n \ 



(109) 



It is instructive to re-derive this solution in a different 
manner, without pre-supposing the form of the solution. 
We look for general self-similar solutions of the form 



p(r,t)=po(t)f(r/r (t)). 



(110) 



We require that all the m ass is in the profile flllO| ) and 
impose the relation (|105| ), implying that 



+oc 



f(x)S d x d - l dx = l 



(111) 



Subs tituting the ansatz (|110|) with Eq . ( P-05D in Eq. 
(103), and imposing the condition (107) where Z is for 
the moment an arbitrary constant, we obtain the differ- 
ential equation 



1 d ( x d-lp-itf_ 

x d ~ 1 dx \ dx 



-2Z l -<{xf + df). (112) 



Noting the identity x d 1 (xf'+df) = (x d f)', this equation 
can be integrated into 



dx 



(113) 



This first order differential equation can again be readily 
integrated. We can choose the constant of integration so 
as to obtain a solution of the form 



/(*)=|[i-(7-iy]i /(7 - l) - 



(114) 



Finally, the normalization condition (111) implies that Z 
is given by Eq. ( |106| ). It is interesting to realize that the 

q-cxponential function e q (x) = [1 + (q — l)x]/ in- 

troduced in the context of Tsallis generalized thermody- 
namics stems from the simple differential equa tion ( |113| ) 
related to the anomalous diffusion equation (103). In- 
deed, the scaling solution of this equation can be written 



f( x ) = -7^ e i{- x2 )i 



(115) 



which generalizes the gaussian distribution obtained for 
the ordinary diffusion equation re cover ed for 7=1. The 
moments (r k ) of the distribution (110) are given by 



(r k )(t)=r (ty 



+ OO 



f{x)x 



fe+d-l 



Sddx. 



(116) 



They exist provided that k > —d for 7 > 1 and provided 
that —d < k < 2/(1 — 7) — d for 7 < 1. They scale like 

{r k ) ocr^J oci«fe/(d+2«). 

The Tsallis entropy is finite for 7 > 73/5 = d/(d + 2) 
and it scales like 



S(t)-nM 



l/n 

-np Q > 



+ OO 



fixySdX^dxotr' 1 !^ 2 ^. 
(117) 



On the other hand, for d > 2, the potential energy W 
-l/(2S d )f(V$) 2 dr scales like 



W ex 



M(r) 

r d—l 



„d-l 



dr 



„rf-2 



,-n(d-2)/(d+2n) 

(118) 
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Comparing Eqs. (117) and (118), we see that the po- 
tential energy is always negligible with respect to the 
entropy for n > n.3. Therefore, the Tsallis free energy 
behaves like 



F(t)+nKM <xt- d/{d+2n \ 



(119) 



for t — ► +00. Note that for 77,3 < n < +00, the free energy 
tends to a finite value —nKM as the system spreads to 
infinity. Alternatively, for the isothermal case n = +00, 
the free energy is given by Eq. (95) of [^6| and it tends 

to — 00. 



We can use the identity (B4) to derive the first cor- 
rection in the evolution of the moments (r fc ) due to self- 
gravity To th at purpose, we introduce the zeroth order 
solution (|104[) in the equation 



d(r k ) 
dt 



k(k + d-2) J Pr k - 2 dr 



r+00 
-fc / r k ~ 
Jo 



dM 
M(r)^-dr. (120) 

or 



The first term gives, after integration, the pure anoma- 
lous scaling 



(r k ) oc t nk/{d+2n) 



(121) 



The second term gives, after integration, the first correc- 
tion due to gravity. If we write A(r k ) = (r ) — {r }q, we 

get 



A(r k ) ex t^T2^+\ 



(122) 



Let us consider some particular cases: (i) for n — > +00, 
we obtain A{r k ) oc ^( fc - d )/ 2 + 1 . If we furthermore con- 
sider the second moment k = 2 (moment of inertia), we 
recover the scaling A(r 2 ) ex t 2 ~ d > 2 of [fUJI. (ii) for k = d, 
we find that A(r d ) = —(d/2)t ex t whatever the index n 
and the dimension of space d. (iii) For n = 77,3, gravi- 
tational effects are of the same order as diffusive effects 
and (r k )o ex A(7 ,fc ) ex t k l d . This case will be studied in 
detail in the next section, (iv) Finally, let us introduce 
the number fco = d — d/n — 2. For k < ho, A(r k ) — ► 0; 
for k = fco, A{r k ) ex 1/t; for k > k , A{r k ) — ► +00. 



B. The critical case n 



II:', 



Finally, for n = 77,3, and since a critical C exists irre- 
spectively of the presence of a confining box, the system 
collapses for O < O c and evaporates for 9 > O c . In the 
latter regime, gravity remains relevant and evaporation 
is controlled by both gravity and the diffusion process 
(see Fig. [15]). Mathematically, this arises from the fact 
that there is an evaporation scaling solution for which 
all the terms of Eq. ( |65|) scale in the same way. More 
specifically, we expect an evaporation density profile of 
the form 



p(r,t) = po(t)f 



ro(t) 



Po(t)r (t) d = 1. (123) 



The relation between the evaporation length r$ and po 
(proportional to the central density) is obtained by re- 
quiring that the diffusive term and the drift term in Eq. 
( p5| ) scale the same. In terms of the mass profile, we have 



M(r,t)=g 



ro(t) 



ith g(x) = / f(x')SdX 
Jo 



ld-1 



dx', 



and in terms of the function s, we have 



(124) 

(125) 

We require that all the mass is contained in the self- 
similar profile teOj, which implies that 



s(r,t) = p (t)S 



ro(t) 



with 



S(x) 



g(+oo) 



+ X 



f{x)S d x d - 1 dx = 1 



(126) 



Inserting the ansatz (125) in Eq. (|70|), using Eq. ( 123 ) 
and imposing 



1 dp 
p% dt 



-dQ, 



d-l dr 

dt 



e, 



(127) 



we obtain the sca ling equation (note the change of sign 
compared to Eq. (\n\)) Q: 



s" + ^±As' + (xS' + ds) 2 ' d fis 

x V 6 



1=0. (128) 



The evaporation radius is given by 

r (t) = (dQt) l ' d . 

The moments scale like (r k ) ex r§ ex 

free energy scales like F(t) + n 3 KM ex i~( d_2 )/ d . 

If we consider the large temperature limit 0»1 where 
the diffusion term dominates on the gravitational drift, 
the foregoing differential equation reduces to 



(129) 
(d&t) k / d and the 



S" 



-S' + (xS' + dS) 2/d = 0. 



In terms of the function / it can be written 

r 2/d f + 



s. 



(d-2)/d 



(130) 



(131) 



which is consistent wi th E q. ( |1 1 3[) u p to the changes 
of notations in Eqs. (107) and ( |l27| ). We can either 
solve t his e quation and impose the normalization con- 
dition (126) or make simple trans form ations in order 
to directly use the results of Sec. VIA. Indeed, let us 
set po = op* and tq = /jr*. We impose p*r d = 1 
leading to ap d = 1. On the other hand, we choose 



a = 2{Sd/Z) ( - d ~ 2 '>/ d where Z is defined by Eq. ( p06| ) so 
that p, = -2d(S d /Z^ d ~ 2 ^ d epl Then, p = ^/,(r/r») 
with /*(x) = crf(x/p). Now, p*, r* and /» have been 
defined so as to coincide with the functions po, ro and / 
of Sec. VIA. Thus, we get f(x) = (l/c)/*(/xx) where /* 
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is the f uncti on (114). Therefore, the normalized solution 
of Eq. (131) with the present notations can be written 



with 



™ = ~Z 



<7fl 



d-1 



2 2 
-fi X 



On 



d/(d-2) 



(d-2)/d 



(132) 



(133) 



and where Z is given by Eq. (106). Proceeding alon g the 
lines of Q, we could expand the solutions of Eq. (12S ) 
(or of the equivalent equation for /) in powers of 8 Mn 
the limit O — > +oo. 




FIG. 15: In d = 3 and for n — ns = 3, we compare the 
evaporation profiles at different times for O = 1 > C , for 
self-gravitating particles (full lines), to the faster evaporation 
dynamics when gravity is switched off (dashed lines). 

In Fig. |l6|, we show the form of the evaporation density 
profile / as a function of > C . As approaches 9 C , 
the central density diverges, whereas the profile tends to 
the one corresponding to free diffusion for large 0. In 
addition, we present numerical simulations for an inter- 
mediate 0, showing that dynamical scaling is perfectly 
obeyed. Moreover, when — > C , we find that the scal- 
ing function obeys itself a scaling relation (see insert of 
Fig. |l|). Defining e = (0 - C )/0 C , wc find 



f(Q,x) 



- l F{x/e 



i/d\ 



(134) 



where F takes the form of a steady polytropic profile of 
index n%. This scaling relation implies that close to C , 
the d-th moment of r scales as 



(r rf (i))~(0-0 c )i, 



(135) 



which is a generalization of our exact result for d = 2 

(7i 3 = +co, t c = i/4) pa, 
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FIG. 16: In d — 3 and for n — ns = 3, we compare the 
scaling profiles for 6 = 0.21 near C « 0.20872, 6 = 1, and 
6 = 100 (top to bottom full lines; for clarity, the 6 = 0.21 
profile has been scaled down by a factor 150). For 9>1, the 
invariant profile corresponds to the Barenblatt solution (pure 
anomalous diffusion) which is a Tsallis distribution with index 
74/3 = 1 + 1/^3. For 6 — > C the invariant profile tends 
to the profile of a steady polytrope with index n^. For an 
intermediate temperature = 1, we illustrate the perfect 
observed data collapse by plotting rg(t)p(r, t) as a function 
of r/r (t), for t = 1.5™ (n = 0, ..., 13). These 14 curves are 
indistinguishable from the theoretical scaling profi le. In the 
insert, we illustrate the scaling relation of Eq. (134) obtained 
for different values of e = (0 — C )/0 C — + 0. 



VII. ANALOGY BETWEEN THE LIMITING 

MASS OF WHITE DWARF STARS AND THE 

CRITICAL MASS OF BACTERIAL 

POPULATIONS 

The generalized Smoluchowski-Poisson (GSP) system 
describing the dynamics of self-gravitating Langevin par- 
ticles shares many analogies with the generalized Kellcr- 
Segel (GKS) model describing the chemotaxis of bacterial 
populations. Below, we briefly review the basic equations 
of chemotaxis and show the close link with the present 
work. 



A. The generalized Keller-Segel model 

The original Keller-Segel model has the form pC 
dp 



c)t 



V • (D 2 (p, c)Vp) - V • (D^p, c)Vc), (137) 



e-^ = -k(c)c + h(c)p + D c Ac. 
at 



(138) 



<r 2 (i))=4(T-T c )t + (r 2 (0)>. 



(136) 



The drift-diffusion equation (137) governs the evolution 
of the dens ity o f bacteria p(r, t) and the reaction-diffusion 
equation ( jl3§| ) governs the evolution of the secreted 
chemical c(r, i). The bacteria diffuse with a diffusion 
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coefficient D 2 and they also move in a direction of a pos- 
itive gradient of the chemical (chemotactic drift). The 
coefficient D\ is a measure of the strength of the influ- 
ence of the chemical gradient on the flow of bacteria. On 
the other hand, the chemical is produced by the bacteria 
with a rate h(c) and is degraded with a rate k(c). It also 
diffuses with a diffusion coefficient D c . In the primitive 
Kellcr-Scgel model, D\ = D\(p,c) and D 2 = D2(p,c) 
can both depend on the concentration of the bacteria 
and of the chemical. This can take into account micro- 
scopic constraints, like close-packing effects pB, U% p3] 
or anomalous diffusion [pi). 

If we assume a constant diffusion coefficient D 2 = D 
and a constant mobility D\/ p = \ (we also consider a 
constant production rate A and a constant degradation 
rate k 2 of the chemical) , we obtain the standard Kcllcr- 
Segel (KS) model 



dp 

at 



V • (DVp - xpVc) 



dc 
l dt 



Ac- k 2 c + Xp. 



(139) 



(140) 



If we now assume that the diffusion coefficient and the 
mobility depend on the concentration of the bacteria, and 
if we set D 2 — Dh(p) and D\ = X9(p)i where h and g are 
positive functions, we obtain the generalized Kcller-Segcl 
(GKS) model M f| F|: 



dp 

at 



^w-(Dh(p)Wp- X g(p)^c), 



9c , 9 

e— = Ac- k 2 c + Xp. 
ot 



(141) 



(142) 



Equation (141) can be viewed as a nonlinear mean field 
Fokker-Planck (NFP) equation BQ] associated with a 
stochastic process of the form 



is a generalized entropy. A straightforward calculation 
shows that 



F 



xg(p) 



= -y f(-Ac + k 2 c-Xp) 2 dr 
{Dh(p)Vp - x.9(p)Vc) 2 dv < 0, (146) 



which is the expression of the i7-theorem in the canonical 
ensemble adapted to dissipative systems. If we consider 
the particular case of a constant mobility g{p) = p and a 
power law diffusion h(p) = 7 / o 7 ~ 1 , with 7 = 1 + 1/n, we 
obtain the polytropic Keller-Segel model |5l||: 



^ = V • (DVp< - xpVc) . 



dc , o 

e— = Ac- k z c + Xp. 
ot 



(147) 



(148) 



The standard Keller-Segel model is recovered for 7 = 
1. Finally, if we neglect the degradation of the chemical 
(k = 0) and consider a limit of large diffusivity of the 
chemical (implying e = 0), we obtain for sufficiently large 
concentrations (see Appendix C of [[73|): 



dp 

7)t 



V • {DVp^ - xpVc) . 



Ac 



-Xp. 



(149) 



(150) 



These equations are isomorphic to the generalized 
Smoluchowski-Poisson (GSP) system (p3|)-(p4|) provided 
that we set 

D = K/t x=Vfc c = -$, X = S d G. (151) 

Therefore, the results of the present paper apply to the 
chemotactic problem provided that the parameters are 
properly re-interpreted. 



^=x(p)Vc +v ^D(^R(i), (143) 

with a diffusion coefficient D{p) = (D/p) J p h{p')dp' and 
a mobility x{p) = X9{p)l P- These equations are asso- 
ciated with a notion of effective generalized thermody- 
namics i49|, |50f| - The Lyapunov functional of the NFP 



equation (141)- (142) can be written in the form of a gen- 
eralized free energy F = E — T e ffS where 



E = Txj^ 



k 2 c 2 ] 



dv 



pedr, (144) 



is the energy, T e ff = D/x is an effective temperature 
given by an Einstein-like relation and 



,5' 



C(p)dr, C"(p) 



Hp) 

9(PV 



(145) 



Formulation of the results with the biological 
variables 



In the gravitational context, we usually fix the coeffi- 
cients £, G and M and use the temperature as a control 
parameter. In the biological context, the coefficients D, 
X and A are assumed given and the control parameter is 
the mass M. Therefore, it may be useful to briefly refor- 
mulate the previous results in terms of the mass, using 
notations adapted to the chemotactic problem. 

For the critical index n = 713 = d/(d — 2) in d > 2 , the 
steady states (polytropes) of the GKS model (149)-(15C) 
exist, in an unbounded domain, for a unique value of the 
mass given by ]59J| : 



M c = S d 



D(l + n 3 ) 



XA 



"3/ ("3-1) 



(152) 
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For d = 3, we have 

M c = 32ttlo 3 I — J -202.8956... ( — J . (153) 

For d = 2, using the identity (|5(j), we recover the critical 
mass 



M n 



8ttD 



(154) 



associated with the two-dimensional standard Keller- 
Segel (KS) model (see |6C]] and references therein). It 
is convenient to introduce rescaled variables so that 
D = A = x = 1- With this system of units the criti- 
cal mass is M c (d) = S d (l + n 3 ) n3/(n3 - 1) w„ 3 = S d [2(d - 
l)/(d-2)] d / 2 u d/{d _ 2) . For example, M c (d = 2) = 8?r and 
M c (d = 3) = 32ttlu 3 = 202.8956.... Using the approxi- 
mate expression of w„ obtained in Eq. (B72) of [plj , we 
can derive an approximate expression of the critical mass 
in the form 

M approx, d j = ^[rf( d + 2)] d ' 2 . (155) 

a 

For d = 2, it returns the exact result M* pprox {2) = M c = 
8tt. On the other hand, M^ rox {d = 3) = 243 and 
M approx( d = 4) = 2842. Using S" d = 27r d / 2 /T(d/2) we 
find that M^ rox {d) ~ 27r d / 2 d d /r(d/2) for d -► +oo. 

Let us briefly discuss the critical dynamics of the GKS 
system with index n = n 3 = d/{d — 2) for d > 2, de- 
pending on the total mass of the bacteria. For M < M c , 
a box-confined system tends to an incomplete polytropc 
confined by the walls of the box. In an unbounded do- 
main, the syst em ev aporates in a self-similar way as dis- 
cussed in Sec. 



VI B 



For M > M c , the system undergoes 
finite time collapse as discussed in Sec. M In a finite time 
t = t co a, it forms a Dirac peak containing a mass M c sur- 
rounded by a collapsing halo evolving quasi self-similarly 
with a time-dependent exponent a(t) tending extremely 
slowly to a — d as t — > t co u. Thus, 



p(r,t)^M c S(r)+ X (r,t), 



(156) 



where x( r ) behaves roughly as r~ d for r — > 0. For 
M = M c , the situation is delicate and depends on the 
dimension of space. For d = 2, in a bounded domain, the 
steady state of the KS model is a Dirac peak (p = +oo). 
We have constructed in jl^] a self-similar solution tending 
to this Dirac peak in infinite time. The central density 
increases exponentially rapidly. In an infinite domain, 
the KS model admits an infinite family of steady state 
solutions parameterized by their central density but the 
Dirac peak (po = +oo) is selected dynamically (the other 
solutions have an infinite moment of inertia and, since 
the moment of inertia is conserved when M = M c , they 
cannot be reached from an initial condition with a fi- 
nite moment of inertia). We have constructed in 16 a 



self-similar solution tending to this Dirac peak in infinite 
time (and ejecting a small amount of mass at large dis- 
tances so as to satisfy the moment of inertia constraint). 



The central density increases logarithmically rapidly. For 
d > 2 and AI = M c , in a bounded domain, the GKS 
model admits an infinite family of steady state solutions 
parameterized by their central density or, equivalently, 
by their natural radius i?*. We have found numerically 
that the system tends to the polytrope where the density 
reaches zero at the box radius (i?* = R). 

Due to the analogy between gravity and chemotaxis 
||53[ , we find that the critical mass of bacterial popula- 
tions in the standard Keller-Segel model in d = 2 and 
in the generalized Keller-Segel model in d > 2 for the 
critical index n = n 3 shares some resemblance with the 
Chandrasckhar mass of white dwarf stars. For example, 
the curves of Figs. and || also represent the mass of 
the bacterial aggregate as a function of the central den- 
sity. As we have seen, they are strikingly similar to the 
mass-central density relation of white dwarf stars in Fig. 
0. Therefore, bacteria and white dwarf stars share deep 
analogies despite their very different physical nature [p9[ . 



VIII. CONCLUSIONS AND PERSPECTIVES: 

THE GSP SYSTEM WITH A RELATIVISTIC 

EQUATION OF STATE 

In this paper, we have studied the critical dynamics, at 
the index n — n 3 , of the GSP system and GKS model de- 
scribing self-gravitating Langevin particles and bacterial 
populations. This study completes our previous investi- 
gation |SlJ that was restricted to the cases n < n 3 and 
n > 7^3. We have seen that, at the index n = n 3 , there 
exists a critical mass M c (independent on the size of the 
system) that is connected to the Chandrasckhar limiting 
mass of white dwarf stars [|59| . In order to strengthen this 
analogy it would be interesting to study the GSP system 
(O)-(H) with the equation of state (0) corresponding to 
relativistic white dwarf stars. In fact, we can already de- 
scribe qualitatively the behavior of the solutions by using 
the results obtained here for polytropes (see also the sta- 
bility results obtained in |6l| ] for relativistic white dwarf 
stars). 

For d = 1 and d = 2, there exists an equilibrium 
state (global minimum of free energy) for all values of 
the mass M. Therefore, the GSP system relaxes towards 
that steady state. 

For d — 3, there exists a critical mass Mchandra = 
0.196701. ..{hc/Gf/ 2 /{iiH) 2 . For M < Mchandra, the 
GSP system tends to a partially relativistic white dwarf 
star (global minimum of free energy). For AI <C 
Mchandra, the density is small so that the equation of 
state reduces to that of a polytropc of index n = 3/2 
(classical limit). Therefore, the GSP system relaxes to- 
wards a classical white dwarf star as described in Fig. 21 
of pl| . For M = Mchandra the density becomes large 
so that the equation of state reduces to that of a critical 
polytropc of index n = 3 (ultra-rclativistic limit). We 
expect that the GSP system forms a Dirac peak of mass 
Mchandra in infinite time. For M > Mchandra, there is 
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no equilibrium state and the system collapses. When the 
density reaches high values, the system becomes equiva- 
lent to a polytrope of index n = 3. Therefore, according 
to the present study, it forms in a finite time a Dirac 
peak of mass Mchandra surrounded by a halo evolving 
quasi self-similar ly with an exponent a(t) converging very 
slowly to a = 3. 

For d = 4, there exists a critical mass M c = 
0.0143958. ..h 4 /{m 2 G 2 p 3 H 3 ) discovered in @. For M < 
M c , the steady states are unstable and the system can 
cither collapse or evaporate (depending on the form of 
the initial density profile and on the basin of attraction 
of the solution) . In case of evaporation, when the density 
reaches low values, the system becomes equivalent to a 
polytrope of critical index n 3 / 2 = n^ = 2 (classical limit). 
In that case, it undergoes a sel f-simi lar evaporation sim- 
ilar to that described in Sec. VI B where diffusion and 



gravity scale the same way. In case of collapse, when the 
density reaches high values, the system becomes equiv- 



alent to a polytrope of index n' 3 = 4 > 



"3 



2 (ultra- 



relativistic limit). In that case, it undergoes a self-similar 
collapse similar to that described in |51j ]. For M > M c , 
there is no steady state and the system collapses in the 
way discussed previously (energy considerations devel- 
oped in [pTJ show that there is no evaporation in that 
case). 

For d > 5, there is no steady state and the system can 
either collapse or evaporate. In case of evaporation, when 
the density reaches low values, the system becomes equiv- 
alent to a polytrope of index n 3 /2 > n^. In that case, it 
undergoes a self -simi lar evaporation similar to that de- 
scribed in Sec. VIA where gravity becomes asymptot- 
ically negligible. In case of collapse, when the density 
reaches high values, the system becomes equivalent to a 
polytrope of index n 3 > 713 . In that case, it undergoes a 
self-similar collapse similar to that described in |5l]] . 

As we have already mentioned, the real dynamics of 
white dwarf stars is not described by the GSP system, 
but is much more complicated. However, we think that 
the study of this simple dynamical model is an interesting 
first step before considering more complicated models. 
At least, it reveals the great richness of the problem. A 
next step would be to take into account inertial effects 
and study the (generalized) Kramcrs-Poisson system and 
the corresponding hydrodynamic equations |17[ . 



Wa = -GM 2 /2 and for d ^ 2, Wu = (d - 2)W where 
W = (1/2) J pQdr is the potential energy. If the system 
is enclosed within a box, we must add a term —dPbV 
on the right hand side, where Pb is the pressure against 
the box. In the following, we assume that the system is 
unbounded so that Pb = 0. 

For a polytropic equation of state P = K p 1 , with 7 = 
1 + 1/n, the free energy (|3q) can be written 



2n 
F = —E km +W- nKM. 
a 



(A2) 



Therefore, the Virial theorem can be expressed in the 
form 

1 dl dF dW , k s 
-£37 = — + W lt + dKM. A3 

2 at n n 

For the critical index 71 = 713 = d/(d — 2), we get 

\^k = (d ~ 2)F + Wii ~ (d ~ 2)W + dKM ' (M) 



For d 7^ 2, it reduces to 



yf t =(d-2)F + d KM. 



(A5) 



For a steady state (1 = 0), the Virial theorem implies 



F eq 



d-2 



KM. 



(A6) 



We have seen in Sec. Ill C that spherically symmetric 
steady states of the GSP system with 12 = 713 exist for 
a unique value of the mass M = M c (for fixed K) or a 
unique value of the temperature = Q c (for fixed M) 
and form an infinite family of solutions param eteri zed by 
their central density po. According to Eq. (A6), they 
all have the same free energy, independent on the central 
density po. Therefore, thcrmodynamical arguments do 
not allow to select a particular solution among the whole 
family. 

For d = 2, the critical index 713 — > +00 and the equa- 
tion of state is i soth ermal with K = ksT/m. Then, the 
Virial theorem (A4) becomes [|16|: 



2^dt V h 



(A7) 



APPENDIX A: VIRIAL THEOREM AND FREE 
ENERGY OF CRITICAL POLYTROPES 

The scalar Virial theorem for the GSP system reads 



I dJ_ 



2-Efcin + Wu 



(Al) 



where I = J pr 2 dv is the moment of inertia, Ekm = 
(d/2) J Pdv is the kinetic energy of the microscopic mo- 
tion and Wu = — J pr ■ V<I> dr is the Virial. For d = 2, 



with ksTc = GMm/A. For a steady state (I = 0), 
the Virial theorem implies T = T c or M = M c . It di- 
rectly yields the result that unbounded two-dimensional 
isothermal spheres exist for a unique value of the mass 
or temperature. The spherically symmetric solution is 
given by Eq. (|3l|) reading 



p{r) 



Po 



[1 + (irp /M)r 2 



(A8) 



This family of steady solutions is parameterized by the 
central density po- The corresponding mass profile is 
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given by M(r) = J p(r')2irr' dr and the gravitational po- 
tential can be obtained from the Gauss theorem d$/dr = 
GM(r)/r with the gauge condition <&(r) ~ GMlnr for 

r — > +00. This yields 



For k — 2, we recover the moment of inertia. Taking 
the time derivative of Eq. (Bl), using the generalized 
Smoluchowski equation (O) and integrating by parts, we 
obtain 



M(r) = 



npor 



1 + {irp /M)r 2 ' 



_ , GM , / M 

$(r) = In h r' 

2 V^Po 



(A9) 



(A10) 



From these expressions, we find that the potential energy 

is 



W 



GM 2 



1+ln 



\TtP0 



(All) 



On the other hand, the Boltzmann entropy Sb 
~ks J(p/m) hi(p/m)dr can be written 



S B = 2Nk B 



2 \m) 



(A12) 



Therefore, the Boltzmann free energy Fb — W — TSb is 
given by 



F B = — 



GM 2 



1 + 1 » { N 



(A13) 



We conclude that the free energy of unbounded isother- 
mal spheres in two dimensions is independent on the cen- 
tral density po- 



I dh 

k K dt 



r k - 2 r-WPdr- / r k - 2 p r • V$ dr. (B2) 



Integrating by parts the first term, we get 



!«§-(*+'<-» 



fc-2 



dr 



„fc-2 



pr ■ V$ dr. 

(B3) 

If we take into account boundary effects, we need to in- 
troduce a pressure term — § Pr k ~ 2 r ■ dS on the r.h.s. 
For k = 2, we recover the Virial theorem (Al). On the 
other hand, for a spherically symmetric system, using the 
Gauss theorem, the second integral can be simplified and 
we obtain 



l^ = ( k +d-2) [ Pr k - 2 dr-G t+ 
k dt . 



M(r)—-dr. 
or 

(B4) 

For k — d, the second integral can be calculated explicitly 

and we get 



APPENDIX B: AN EQUATION FOR THE 
MOMENTS (r k ) 



^ = 2(d-l) / Pr d - 2 dr 
dt J 



d 



GM 2 



(B5) 



Let us introduce the moments of order k: 
h{t)= f pr k dr. 



For c2=l, the first term on the r.h.s. must be replaced 
(Bl) by2P(0,t). 
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Since the free energy F[p] coincides with the energy func- 
tional W[p] of a barotropic gas (up to a positive macro- 
scopic kinetic term) J65[ , we conclude that p e <j(r) is lin- 
early dynamically stable with respect to the GSP system 
iff it is formally nonlinearly dynamically stable with re- 
spect to the barotropic Euler-Poisson system p5|. 
We shall study the problem in d dimensions because: (i) 
We have found that the structure of the mathematical 

troblem with the dimension of space is very rich WA M, 
y], exhibiting several characteristic dimensions, (ii) In 
grav ity, t he usual dimension is d = 3 but, in biology (see 
Sec. VII), the bacteria (or cells) are compelled to lie on 
a plane so that d — 2. 

It may appear artificial to put the system in a "box" . In 
gravity, the box delimitates the region of space where the 
system can be assumed isolated from the surrounding and 
whe re statistical mechanics applies. In biology (see Sec. 



VII), the box has a physical meaning since it represents 
the container in which the bacteria (or cells) are confined. 
The reader should be aware that, in the sections dealing 
with the dynamics, po and ro do not exactly coincide 
with the quantities of the same name introduced in the 
sections dealing with the statics (they usually differ by a 
factor of proportionality). 

More generally, for a polytrope of index n we have p oc 
r~ a at t — tcoii, with a = 2n/(n — 1) so that the self- 
similar solution exists provided that a — d+1 < 1 leading 
to 1/n < 1/n-s (i.e. n > nz for d > 2). This is precisely 
the range of indices for which the complete polytropes 
are dynamically unstable [ pl[ . 

Defining p c (x) as the equilibrium density profile at O = 
<d c , then the first component can be written pi(r,t) = 
^Pc(r/r ). 

Looking for a self-similar solution of the form ( |7l| ) for 
any index n and requiring that the diffusion and gravity 
scale the same, we find that po?"o ~ 1 where a is given by 
Eq. ([72). The profile will contain all the mass provided 
that po^o ~ 1 which implies a — d leading ton = 713. 
Thus, it is only for the critical index that we can have a 
self-similar solution where the diffusion and gravity scale 
the same and which contains al l th e mass. 
The scaling equations ( |77| ) and ( 12S ) have a very different 
mathematical structure. The scaling equation for collapse 
(|77|), valid for n > 713, leads to an eigenvalue problem for 
S(x) I12I pi I . Indeed, it admits a physical solution for 
a unique va ue of S(0) equal to S, (say). For 5*(0) < 
S*, the solution becomes negative at some point, and for 
S(0) > 5**, it diverges at a finit e 3*0 . By contrast, the 
scaling equation for evaporation (125), valid for n = 713, 



admits a one parameter family of solutions parameterized 
by 5*(0). Then, the suita ble value S* is selected by the 
normalization condition ( 12£ ). 



